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' The traffic of molecular motors tfirougfi open tube-like compartments is studied using lattice 

, models. These models exhibit boundary-induced phase transitions related to those of the asymmetric 

fSJ . simple exclusion process (ASEP) in one dimension. The location of the transition lines depends on 

» I ' the boundary conditions at the two ends of the tubes. Three types of boundary conditions are 

CIh, studied: (A) Periodic boundary conditions which correspond to a closed torus-like tube. (B) Fixed 

• motor densities at the two tube ends where radial equilibrium holds locally; and (C) Diffusive motor 

injection at one end and diffusive motor extraction at the other end. In addition to the phase 
diagrams, we also determine the profiles for the bound and unbound motor densities using mean 
field approximations and Monte Carlo simulations. Our theoretical predictions are accessible to 
experiments. 
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I. INTRODUCTION 



Molecular motors are proteins that transform the free energy released from chemical reactions into mechanical 
work. In this article we consider a special class of motor proteins, namely cytoskeletal motors which perform directed 
^ walks along cytoskeletal filaments as reviewed in [1,2]. In the cell, these motors have different functions related, e.g., 
to vesicle transport and cell division. The best studied examples are kinesins, which walk along microtubules, and 
certain types of myosins, which walk along actin filaments. After a certain walking time, such a motor unbinds from 
C ■ its filament because its binding energy is finite and can be overcome by thermal activation. For kinesins, this typically 
] happens after 100-150 steps or after a walking time of about 1.2-1.8 seconds, see, e.g., [3,4]. In many motility assays, 
(~| the filaments are immobilized on a substrate and are in contact with an aqueous solution. In such a situation, the 
O unbound motors diffuse in the surrounding fluid until they eventually reattach to the same or another filament. 
, O , . Recently, we have introduced lattice models to study the motors' random walks, which consist of many diffusional 
encounters between motors and filaments in open and closed compartments [5,6]. If many motors are placed in such 
a compartment, hard core exclusion between the motors has to be taken into account, since the motors are strongly 
attracted to the binding sites of the filaments, so that the filaments get overcrowded. These models are new variants of 
driven lattice gas models and exclusion processes, for which the active processes which drive the particles are localized 
to the filaments. 

Lattice models of driven diffusive systems have been studied extensively in the last years, see e.g. [7,8]. The simplest 
: model is the asymmetric simple exclusion process (ASEP) in one dimension, where particles hop on a one-dimensional 
lattice with a strong bias towards one direction (in the simplest case there are no backward steps at all) and the only 
interaction of the particles is hard core exclusion, i.e., steps to occupied lattice sites are forbidden. When coupled to 
^ open boundaries, this simple model already exhibits a complex phase diagram, see e.g. [8], which we will review below 
in some detail. The first model for the 1-dimensional ASEP was introduced more than 30 years ago by MacDonald 
^ • et al. [9,10] in the context of protein synthesis by ribosomes on messenger RNA (mRNA). At that time it was solved 
using a mean field approach and used to explain results of radioactive labeling experiments [11-13] which showed 
that protein synthesis gets slower as the ribosome moves on the mRNA template. The model of MacDonald et al. 
explained this by the steric hindrance between successive ribosomes along the mRNA track. Two years later the 
Q , same model was discussed by Spitzer as a simple example for interacting particles in probability theory [14]. Since 
then, the asymmetric simple exclusion process (and variants) has been studied extensively as a generic model for 
non-equilibrium phase transitions [15-17] and interacting stochastic systems [18] as well as in other applications such 
as traffic fiow [19]. Many properties of the 1-dimensional ASEP are known exactly. 
5h ' As mentioned, the lattice models for random walks of molecular motors differs from the driven lattice gas models 
in an important way: Walks of molecular motors are only 'driven' as long as the motor is attached to a cytoskeletal 
filament, hence 'driving' is localized to one or several lines. It can therefore be viewed as an ASEP which has the 
additional property, that particles (molecular motors) can unbind from the track with a small probability, diffuse in 
the surroundings and reattach to the same or another filament. In more mathematical terms, the ASEP is coupled 
to a symmetric exclusion process via adsorption and desorption of particles onto filaments. 

Boundary conditions play an important role in driven systems. This becomes apparent, e.g., if one compares a 
tube-like system with periodic boundary conditions with one with closed boundaries. In the system with closed 
boundaries, a traffic jam of motors arises at one end of the system and the current of motors bound to the filament 
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is balanced by diffusive currents of unbound motors as first shown in [5]. With periodic boundary conditions, motors 
arriving at the right end of the system just restart their walk from the left end and a net current through the systems 
is obtained. 

In this article, we study the stationary states of tube-like compartments with open boundaries. These compartments 
have the shape of a cylinder and contain one filament which is placed along the cylinder axis in order to obtain the 
simplest possible geometry. The bound motors move along the filament and the unbound motors diffuse within the 
cylinder, see Fig. 1. At the ends of the tube, motors are inserted and extracted. Such a system is accessible to in vitro 
experiments using standard motility assays, but it can also be viewed as a strongly simplified model for motor-based 
transport in an axon [20], if these motors, which are synthesized in the cell body, are at least partly degraded at the 
axon terminal [21], a situation that can be mimicked by insertion and extraction of motors at the ends of a tube. The 
stationary states depend strongly on the way, in which the motors are inserted and extracted at the boundaries as we 
will explicitly demonstrate for three different types of boundary conditions, see Fig. 2. 

Our article is organized as follows. After introducing the model in section II, we start in section III with periodic 
boundary conditions. This case can be solved exactly, since it satisfies local balance of currents in the radial direction, 
see appendix A. In section IV, we discuss the situation in which the density of bound motors on the filament is 
fixed at the boundaries. Finally, in section V, we consider the case where the filament is shorter than the tube and 
the motors diffuse into and out of the tube. The main tool to study the open systems are Monte Carlo simulations. 
These are supplemented by dynamical considerations and self-consistent or mean field calculations. Some details of 
the latter calculations are presented in appendix B. 

II. THEORETICAL MODELLING 
A. Tube geometry 

We consider the motion of molecular motors in a cylindrical tube as shown in Fig. 1. The tube has length L and 
radius R. The total number of motors within the tube, denoted by A^mo, defines the overall motor concentration 

Note that this concentration corresponds to the particle number density of the motors and, thus, is insensitive to the 
size of the motor particles. In general, the latter size depends on the type of motor and on the type of cargo attached 
to it. In the following, we will implicitly assume that the motor particle has a linear size which is comparable to 
the basic length scale £ as defined further below. If one wants to study the dependence on the motor particle size in 
a systematic way, one should measure the overall motor concentration in terms of the volume fraction of the motor 
particles. 

The cylindrical tube contains one filament located along its symmetry axis which is taken to be the x-axis; the 
two other Cartesian coordinates are denoted by y and z. Motors bound to the filament undergo directed motion, 
while unbound motors diffuse freely. Since the motors are strongly attracted by the filament, a large fraction of these 
motors is in the bound state and mutual exclusion from the binding sites on the filament has to be taken into account 
even for relatively small overall motor concentrations. 

In order to include this mutual exclusion (or hard core repulsion) in the theoretical description, we map the system 
onto a lattice gas model on a simple cubic lattice. The lattice is oriented in such a way that its three primitive vectors 
point parallel to the x-, y-, and 2-axis, respectively. A rather natural choice for the lattice parameter i is the repeat 
distance of the filament, which is 8nm in the case of kinesin motors moving on microtubule and 36 nm for myosin V 
motors moving on actin filaments. 

The discretized tube consists of one line of binding sites, which represents the filament, and Nch unbound 'channels', 
i.e. lines of lattice sites parallel to the filament. Thus the cross section ^ of the tube is equal to 

</>=(l+AehK' . (2) 

For sufficiently large radii, one has </> « 7ri?^, while for small radii, there are corrections due to the underlying lattice. 

In the following, we will measure all distances such as L and R in units of £. Thus, for the simulations, both L 
and R will be quoted as integers. The integer value of R corresponds to a certain number of lattice sites along the 
Cartesian coordinates y and z which are perpendicular to x and run parallel to two basis vectors of the simple cubic 
lattice. Thus, the interior of the tube contains all channels with y'^ + <R?. 
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B. Random walks with mutual exclusion 



At each time step, a bound motor attempts to make a forward or backward step and to jump to the next lattice site 

to its right or to its left with probability a or /?. In addition, the motor attempts to jump to each of the neighboring 
sites away from the filament with probability e/6, and does not attempt to jump at all, i.e., to rest at the filament 
site, with probability 7. Since the sum of these probabilities is equal to one, we have 7 = ! — a — /3 — 2e/3. The 
velocity Wb in the boimd state is given by Wb = (a — /3) ^/t, where r is the basic time scale of these random walks. In 
the following, we will measure all times and rates in units of r and t^^, respectively. Since backward steps are rare 
for cytoskeletal motors, we will focus on the case /3 = 0, i.e., we will ignore backward steps. 

If the particle unbinds from the filament, it attempts to jump to all nearest neighbor sites of the simple cubic 
lattice with equal probability 1/6. By choosing the time scale r as i?^/Z)ubj this hopping probability can be made to 
fit the diffusion coefficient of unbound motors, -D,ib. When measiirecl in units of I and r, the dimensionless diffusion 
coefficient i^ub = 1/6- In principle, the resting probability 7 can then be used to account for the ratio D^\^l(v\,l), 
which is quite large in the case of kinesin, and to adapt the velocity v\, to the values obtained from experiments [5] . 
For simplicity, we will often choose 7 = in order to eliminate one parameter from the problem. If an unbound motor 
attempts to hop to a filament site, the motor binds to it with sticking probability TTadj while the step (and hence 
binding) is rejected with probability 1 — TTad- 

Both in the bound and in the unbound state, hopping attempts can only be successful, if the target site is not 
occupied by another motor; otherwise the particle must stay where it is. In the following, we will mainly study overall 
motor concentrations Pmo in the range < pmo ^ 0.05. For such small values of pmo, miitual exclusion of the unbound 
motors can be safely ignored. However, because the motors are strongly attracted to the filament, the concentration 
of bound motors is much larger and it is crucial to take their mutual exclusion into account. 



III. PERIODIC BOUNDARY CONDITIONS 



First, we consider a cylindrical tube with periodic boundary conditions in the longitudinal direction. Because of 
the translational invariance in the direction parallel to the filament, there are no net radial currents in the stationary 
state. Indeed, a non-zero radial current in a state, which is translationally invariant in the longitudinal direction, 
would lead to net radial transport of motor particles, which is incompatible with the reflecting radial boundaries. This 
means that there is a bound current on the filament, but both currents of motors binding to and unbinding from 
the filament and radial currents of unbound motors are balanced locally. We call this situation radial detailed balance 
or radial equilibrium. If there is only one unbound channel, A'ch — 1, (or iVch equivalent channels) radial equilibrium 
is equivalent to adsorption equilibrium. It is clear, that this will no longer be true, if translational invariance is broken 
by boundaries or blocked sites on the filament. Another important property of systems with periodic boundary 
conditions is that, in this case, the number A^mo of motors in the system is conserved, which does not apply to open 
systems. 

Because of translational invariance along the x-axis, the bound and unbound motor densities, pb and pubi do not 
depend on a;. In addition, it follows from the absence of radial currents that pub is also independent of the radial 
coordinate r = (y^ + 2^)^/^. Hence pb and pub are constant, and radial equilibrium implies the relation 



which leads to 



|pb(l - Pub) = ^Pub(l - Pb) (3) 

e/TTad + (1 - e/7rad)/0ub 

Here and below, the densities pb and pub are local particle number densities which satisfy 

< Pb < 1 and < Pub < 1 • (5) 

In dimensionful units, this corresponds to < pb < l/£^ and < Pub < 1/^^. 

Since the total number A^mo of motors is conserved, we can use the normalization condition, 

Pb + AchPub = (6) 

to obtain a quadratic equation for the densities as a function of the system size and the total number of motors. Since 
one root is always negative, the physically meaningful solution is: 
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(7) 



From this expression for the unbound density, the bound density pi, follows via Eq. (4) and the stationary eurrent is 
given by J = jb = fb/Ob(l — Pb)- The current calculated in this way is shown in Fig. 3 as a function of N^o/L. The 
data points (circles) are the results of Monte Carlo simulations for a system of length L = 200 and radius i? = 25 and 
are in very good agreement with the analytical solution. 

It follows from the analytical solution that the current J = jb vanishes at N^o/L = 1 + A^ch, behaves as 

^« I ^ (8) 

Vb 1 + ^iVeh L ^ ' 

for small N^-^o/L, and has the maximal value max(J/t;b) = 1/4 for pb = 1/2, pub = -^/(l H — —) and N^no/L = 
l/2 + iV,,^/(l + ^). 
Let us add two remarks: 

(i) Eq. (3), as stated here, can be considered as a mean field equation. However, using the quantum Hamiltonian 
representation of the stochastic process, it can be shown to hold exactly. The calculation is simple, but rather technical 
and is therefore presented in Appendix A, which shows that stationary states can be constructed as product measures 
provided the bound and unbound densities satisfy Eq. (3). 

(ii) Note that in contrast to a homogeneously driven lattice gas such as the asymmetric simple exclusion process, 
there is no particle hole symmetry here. Particles attempt to leave the filament with rate e/6 to a neighboring site 
while holes do so with rate 7rad/6, i.e. particles are strongly attracted by the filament, while holes are not. However, 
if one considers only the bound density, the current density relationship J = jb — VbPb{^ — Pb) is invariant under 
the exchange of particles and holes. As we will see in the next section, this can lead to an apparent particle hole 
symmetry for systems with radial equilibrium, because the radial currents vanish and the state of the system can be 
determined by the bound density alone. 



IV. OPEN BOUNDARIES WITH RADIAL EQUILIBRIUM 

Now, let us consider the more interesting case, where the tube is open and the densities at the left and right 
boundary are fixed. To be precise, we consider two different sets of boimdary conditions, (B) and (C) as shown 
schematically in Fig. 2. In this section, we study case (B) while case (C) will be considered in the next section V. 

For case (B), we add two layers of boundary sites at a; = and x = L + 1 with y^ + z"^ < E?. As before, the filament 
is located at y = ^ = 0. We then fix the density on the additional filament sites according to 

Pb{x = 0) = Pb.in and pb{x = L + l) = pb,ex , (9) 

see Fig. 2. Furthermore, the densities on the nonfilament boundary sites are chosen in such a way that radial 
equilibrium as given by (3) holds at both boimdarics. 

These boundary conditions are implemented by the following choice of random walk probabilities. First, we eliminate 
two parameters from the problem, namely the jump probability /3 to make backward steps on the filament and the 
resting probability 7 to make no step at all on the filament. Thus, we take /3 = 7 = throughout this section. 

Next, when we choose a site within the left or right boundary layer during the Monte Carlo sweep, we first draw 
a random number u) which is uniformly distributed over the interval < w < 1. The chosen left or right boundary 
site is taken to be occupied if w < pb,in or w < Pb.ex, respectively. If a boundary site with j/^ + > is occupied 
by a motor particle, this particle attempts to jump into the tube with probability 1/6. If the left boundary site with 
y = z = is occupied, the corresponding particle attempts to jump onto the left end of the filament with probability 
(1 — 2e/3). If the right boundary site with y = z = is occupied, this particle cannot enter since /3 = 0. In addition, 
all particles which jump from a site within the tube, i.e., from a lattice site with 1 < x < L, onto a boundary site, 
are extracted from the tube. 
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A. Phase diagram 



In the limiting case with A^ch = and e = 0, our system becomes equivalent to the 1-dimensional ASEP, for which 
the density profiles and the phase diagram arc known exactly [16,17]. Let us therefore summarize some of the known 
properties of this process, for which we use symbols without the subscript 'b'. 

For the 1-dimensional ASEP, there are three different phases. If the density pin at the left boundary is small and 
satisfies < 1/2 and if the density pcx at the right boundary is not too large with p^x < 1 — Pm, the system is in the 
low density (LD) phase, for which the bulk density is equal to the left boundary density, and the current is given by 
J = wpin(l — Pin)- Bccausc of the particle hole symmetry, an analogous situation holds for pex > 1/2 and pin > 1 — Pox 
[high density (HD) phase]. Now, the bulk density is given by p° = pox and the current is J = wpex(l — Pex)- At 
the line pox = 1 — Pin < 1/2, a discontinuous phase transition takes place, and the bulk density jumps from the left 
boundary value to the right boundary value. Finally, for pin > 1/2 and pcx < 1/2, the bulk density p*^ = 1/2 and the 
current attains its maximal value J = v/4,. Therefore, this phase is called maximal current (MC) phase. The phase 
transition towards the maximal current phase is continuous with a diverging correlation length. 

The formation of three different phases can be understood in terms of the underlying dynamics of domain walls and 
density fluctuations [8] . In the low density and high density phases, the selection of the stationary state is governed by 
domain wall motion. Thus consider a domain wall, which forms between regions of different densities. If its velocity 
is positive, the domain wall travels to the right boundary and the bulk density is equal to the left boundary density 
(low density phase). Likewise, the bulk density is given by the right boundary density, if the domain wall velocity is 
negative (high density phase). At the transition line with pex = 1 — Pin < 1/2 the domain wall velocity is zero, and 
domain walls diffuse through the system. 

The second mechanism is related to density fluctuations. If a small perturbation of the density, corresponding, e.g., 
to some added particles, enters the system at the left boundary, it can move with positive or negative velocity, i.e., 
it can spread into the bulk or is driven back towards the boundary. In the maximal current phase the velocity of 
such density perturbations, is negative and density fluctuations coming from the left boundary are driven back to the 
boundary. Hence, increasing pjn does not increase the bulk density, since additional particles cannot enter the system 
(overfeeding effect). Both the velocities of domain walls and of density fluctuations are governed by the same density 
current relation which is j = vp{l — p) for the 1-dimensional ASEP [8]. 

For the filament in a tube as considered here, the velocities of domain walls and density fiuctuations are similar to 
those for the ASEP in one dimension. This can be understood from the fact that the density-current relationship on 
the filament is the same as for the 1-dimensional ASEP provided one rescales all currents by the factor (1 — 2e/3) 
which arises from the possibility to unbind from the filament. 

Thus, let us first consider the case, where the behavior on the filament is determined by domain walls, i.e. the 
low density and high density phases. Drift motion of domain walls is governed by the domain wall velocity on the 
filament, which is Vs = v{l — Pex — Pin) for the one- dimensional ASEP [8]. In the tube system considere here, the 
domain wall velocity is slowed down compared to this value, because the domain wall of the unbound density must 
follow the bound density domain wall by binding and \mbinding of motors to and from the filament. However, the 
sign of the domain wall velocity is the same for this tube system as for the one-dimensional ASEP, and the domain 
wall velocity changes sign at the same values of the boundary densities. An explicit expression for Vg can be obtained 
from the general expression for the domain wall velocity as given in Ref. [8] by integrating the density over the tube 
cross-section, which leads to 

VbPb.ox(l - Pb.ox) - l'bPb.in(l - Pb.in) „s 

Vg = '■ '■ '■ (10) 

Pb,cx — Pb,in + ^chpub.cx — ^chPub,in 

for the geometry considered here. Remember that the unbound densities are related to the bound densities by radial 
equilibrium at the boundaries. If a domain wall spreads from the left or right boundary into the system, radial 
equilibrium will hold approximately in the bulk because of translational invariance, but, in addition, all the way down 
to the dominating boundary, where we have imposed radial equilibrium via the boundary conditions. Therefore, the 
current and the bulk density are determined by the bound density at the boundaries. 

The drift velocity of fluctuations of the bound density is Vc = I — 2pb [8]. If the behavior on the filament is 
determined by this velocity, i.e. in the case pb,in > 1/2 and pb,cx < 1/2, radial equilibrium will not hold up to the 
boundary and unbound motors can enter the tube. However, after a short distance they will bind to the filament 
and act like an additional particle in the maximal current phase: The density fiuctuation generated by the additional 
particle moves towards the boundary. Therefore in the bulk, the bound density is 1 /2 and the current is Vb /4. 

Summarizing these considerations, we predict to find exactly the same phase diagram as for the 1-dimensional 
ASEP, see Fig. 4. In fact, we have chosen our boundary conditions (B) in such a way that we only have to replace the 
density p of the 1-dimensional ASEP by the bound density pb and the boundary densities pin and Pex by the boundary 
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densities on the filament, pb.in and pb.ox- The unbound density in the bulk is obtained by radial equilibrium from the 
bound density. These cxpeetations are confirmed (i) by a detailed analysis of the discrete mean field equations and 
(ii) by extended Monte Carlo simulations. 

A remarkable feature of the phase diagram is that it seems to exhibit particle hole symmetry, while the dynamics 
does not. The signature of particle hole symmetry in the phase diagram is the symmetry between the high density 
and low density phase. In both these phases, the bulk density is approximately constant and radial equilibrium holds 
approximately in the whole system except close to the left or right boundary. Radial equilibrium holds exactly at that 
boundary which determines the bulk behavior. Therefore radial currents vanish on average and the phase diagram 
is determined by the bound density alone, resulting in a phase diagram which gives the impression of particle hole 
symmetry, although this symmetry is broken. Indeed, even though the substitution of the bound density pb by 1 — pb 
leads to the same bound current j^, it does not lead to the unbound density 1 — Pub- Two stationary states with 
bound densities that are related through particle hole symmetry are characterized by two unbound densities which 
are both smaller than the bound ones and thus break the apparent symmetry. 



B. Density profiles 

To discuss the concentration profiles of the bound and unbound motors, we use continuum mean field equations and 
compare the mean field results with simulations. In addition, we make a two-state approximation, i.e., we consider 
the case of a single unbound channel, so that a motor can be in only two states, bound to the filament or unbound. 
The two-state model is exact for an arbitrary number iVch of equivalent unbound channels, but it can also serve as an 
approximation for the original tube systems: Since the unbound density depends only weakly on the radial coordinate, 
we will consider the approximation in which p^h is taken to be independent of r and the unbound channels are taken 
to be equivalent. This approximation corresponds to a two-state model in which the bound motors are described by 
the density pb{x) and the unbound motors by the density puh{x). The effective diffusion coefficient for the unbound 
motors is given by NchD^y^. 

Using again the mean field approximation, the total current J parallel to the tube axis is now given by 

d d 

J = ^^bPb(l - Ph) - ^bg^Ph - A^ch-Dub^Pub (11) 



and the equality between incoming and outgoing currents at any lattice site leads to 

d_ 
dx 



d 

■i'bPb(l - Pb) - -Db^Pb = 7radPub(l - Pb) - epb(l - Pub), (12) 



dx 

where Db and i?ub are the diffusion coefficients of motors in the bound and unbound state, respectively. In addition, 
we have introduced the rcscaled binding and unbinding rates, ffad = f ""ad and e = |e. 

In order to determine the density profiles far from the boundaries, we first calculate the homogeneous, x-independent 
solutions, p^ and p\, of the two mean field equations (11) and (12). The first equation (11) for the total current J 
then reduces to the current-density relationship 

J = v^pI(1 - pI) (13) 

whereas the second equation (12) becomes 

i?adPSb(l-pS)=f>S(l-pub) (14) 

which implies radial equilibrium for the x-independent solutions. We then decompose the densities according to 

Ph{x) = Pb + "nhix) and Pub(a;) = p^b + iluh{x), (15) 

and expand the mean field equations (11) and (12) in powers of the density deviations r/b and r/ub- The details of this 
expansion are described in Appendix B. 



1. Low density and tiigti density phases 

For the high and low density phases with pjj 7^ 1/2, the expansion of the mean field equations up to first order in 
r?b and rj-ah leads to an exponential approach ~ exp(a;/^) of the density profiles towards the homogeneous solutions 
p° and p°jj, see Appendix B. The corresponding decay length ^ satisfies the cubic equation 
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-Vb{l - 2pg)^3 + (Db + N,hD^^g)e + jVch^^ub '"'^^ ^ ^^^^ C - A^chi'ub^ = (16) 
which is solved numerically. 

For the ASEP in one dimension, one has A^ch = and (16) reduces to -Wb(l - 2p°)^^ + Db^^ = -^^(^ - Co) = 
with 

^o = Dy,/Ml-2pl) . (17) 

Thus, in this limit, mean field theory leads to the correlation length £, = £,0- For iVch > 0, wc choose the unique 
solution of (16) which approaches (17) as A^ch vanishes. This solution behaves as ^ « Co(l + gNchDuh/D^) for small 
A^ch. 

In addition to the mean field calculation, we again used Monte Carlo simulations in order to determine the density 
profiles as shown in Figs. 5 and 6. As predicted by the mean field calculation, the constant bulk densities for the 
bound and unbound states are approached exponentially in the low and high density phases. The corresponding 
decay length ^ is found to be the same for the bound and the unbound density and to diverge as one approaches the 
maximal current phase. 

Some simulation results for the decay length ^ are displayed in Fig. 7. In this case, the radius R of the tube was 
varied while the boundary densities were kept fixed. The latter densities were chosen to be pb,in = 0.38 and pb,cx = 0.6 
which lies within the low density phase but is close to the phase transition line which separates the low from the 
high density phase. One surprising feature of the Monte Carlo data for the decay length ^ is that they exhibit a 
pronounced minimum as a function of tube radius R. 

For comparison, we also display in Fig. 7 the ^-values as obtained from several mean field approximations corre- 
sponding to the dashed, dotted and solid lines. The dashed line is obtained from the solution of equation (16) which 
we derived from the continuous mean field approximation for the two-state model. We also determined this quantity 
using a lattice version of this approximation (dotted line) and a more elaborate mean field approximation (solid line) 
in which we solved the diffusion equation in the cylindrical compartment and matched this solution to the directed 
transport along the filament. 

Inspection of Fig. 7 shows that all three mean field approximations are quite consistent with each other and lead 
to C ~ i?^ as follows from (16). Such an increase of ^ for large R is in fair agreement with the Monte Carlo data 
displayed in Fig. 7. However, in contrast to the Monte Carlo simulations, all three mean field approximations give a 
monotonic increase of ^ with increasing R. 

The largest discrepancy between the mean field results and the Monte Carlo data is found in the limit of small 
R for which one recovers the 1-dimensional ASEP. In this limit, the Monte Carlo data should be quite reliable as 
one concludes from the value obtained for = which is in very good agreement with the exact solution for the 
1-dimensional ASEP as given by [16] 



1 _^ 1 



with ^fe = -— — 1 (18) 

ln[4pfe(l - pk)\ 



for k = in, ex. Thus, we conclude that the decay length ^ does indeed exhibit a nonmonotonic dependence on the 
tube radius R and that this behavior is not reproduced by the mean field approximation. 

Finally, we note that the different behavior of the decay length £ for large and for small R is correlated with a 
qualitatively different behavior of the corresponding density profiles as observed in the Monte Carlo simulations. As 
an example, let us consider the density profiles within the low density phase as in Fig. 5. In this case, the bound and 
unbound densities exhibit plateau regions which are determined by their values Pb.in and Pub.in at the left boundary. 
As one gets closer to the right boundary where the motors can leave the tube, the densities start to deviate from these 
constant values, and these deviations grow exponentially as ~ exp(a;/C). For small R, the corresponding profiles are 
convex upwards for all values of x. For large R, on the other hand, the profile exhibits an inflection point close to the 
right boundary. This inflection point moves towards the interior of the tube as R is further increased. 



2. Maximal current phase 

In the maximal current phase, one has p'^ = 1/2, and one has to consider terms up to second order in the density 
deviations rj^ and rj^h, see Appendix B. One then flnds that rji, satisfles the nonlinear differential equation 

d 
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with 



,_ g(l-Pub)+W°b 
7^ad(l-plJ)+e>{^ 



This differential equation can be solved by separation of variables. The solution behaves as 

r]b ~ tor large x (21) 

VbX 

i.e., the deviation of the bound density from its asymptotic bulk value = 1/2 decays as ^ 1/x within the mean 
field approximation. 

For Nch ~ 0, this becomes identical with the mean field solution for the ASEP in one dimension as discussed in 
[15]. In this latter case, an exact solution is available [16,17], which shows that the density profile decays as ~ l/x^^^, 
i.e., with a different exponent. This is related to the fact that density fluctuations spread superdiffusively in the 1- 
dimensional ASEP [22]. The dispersion ((Aa;)^) of an ensemble of particles in such a system behaves as ((A.x)^) ^ if^l"^ 
for large times t. This superdiffusive spreading of density fluctuations in one dimension can be taken into account 
within a mean field approximation if one considers a scale-dependent diffusion coeflacient as shown in Ref. [15]. In 
the following, we will extend this approach to the tube geometry considered here. 

Thus, we now replace Db in the mean field equation (19) by a scale-dependent diffusion coefficient D\,{x) and 
consider the modified mean field equation 

[I)b(a:) + <?iVch£'ub]^%~-i^br/b • (22) 
A convenient choice for D\,{x) which embodies the correct superdiffusive behavior for the ASEP in one dimension is 



Xs 



I)b(x) = D,e( . (23) 



The left boundary is located at a; = a;o, and Dsc and x^c represent two scale parameters. With this choice, the 
modified mean field equation (22) can again be solved by separation of variables. As a result, we obtain 



r/b(a;) = a 

where we have introduced the abbreviations 



^-Mn(l + Vy/6) 



-1 



(24) 



%(a;o) 

y={x- xo)/xsc , (25) 



a = Dsc/2vbXs 



(26) 



and 



(27) 



The 'initial' value ?7b(a;o) = Pb,in — 5 denotes the density deviation at the left boundary. 

Far from the left boundary, i.e., for large values of y ~ x, the expression (24) leads to the asymptotic behavior 



a 



HVy) 
Vv 



(28) 



Thus, the deviation of the bound density from its asymptotic value now decays as 1/y^ l/\/^ as for the ASEP in 
one dimension. 

In addition, we also obtain a correction term in (28) which depends on 6 = 1 +g NchDuh/Dsc- For large tube radius 
R, one has A^ch ^ and h ~ A^ch ~ R^- Therefore, the correction term becomes large for large R. We can now 
define a crossover length y = y^, which the correction term has the same size as the leading term. This leads to the 
implicit equation 



(29) 
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and, thus, to 



Close to the left boundary, i.e., for small values of y ~ a; — a;o, the expression (24) for the deviation of the bound 
density from its asymptotic value = 1/2 leads to 

Vh{x) « ??b(a;o)[l - ^'?b(^o)y]. (31) 

We may now define a second crossover length or extrapolation length y = y** at which the two terms in (31) cancel. 
This extrapolation length is given by 

y** = 2ab/T]b{xo) = 2a6/(pb,iii - ^) • (32) 

For large tube radius R, the latter length scale grows as y** ~ 6 ~ ii^. In general, the extrapolation length can have 
both signs but, for the maximal current phase, one always has pb,in > 1/2 and, thus, y,, > 0. 

There is also an intermediate range of y-values defined by y** <C y <C 6^ ~ y*. For these y-values, the deviation -q^ 
of the bound density from its asymptotic value as given by (24) simplifies and becomes 

?7b ~ = • (33) 

y V},Xsc y 

Thus, for these intermediate y-values the density deviation decays again as ~ 1/y. 

In summary, the theory described here indicates (i) that the 'initial' value at the left boundary is felt up to an 
extrapolation length y** ~ 6 ~ R^, (ii) that the true asymptotic behavior of the density deviation is obtained for 
y > y* ^ -R^ as follows from (30), and (iii) that the density deviation decays as 1/y on intermediate length scales 
with y** <C y <C y* . 

These conclusions agree with the results of Monte Carlo simulations. In these simulations, the tube length L is 
necessarily finite. This implies that the profiles observed in the simulations may not reach the asymptotic behavior 
present in a tube of infinite length. Since both the crossover length y* and the extrapolation length y„ increase 
quickly with increasing tube radius R, we expect to find the true asymptotic behavior only for sufficiently small 
values of R. This expectation is confirmed by the simulation results. For sufficiently small R, the density deviation 
77b is found to decay as ~ for large x, as shown in Fig. 9 for i? = 3. For sufficiently large values of R, on the 

other hand, the observed profiles decay as 1/x, see the Monte Carlo data in Fig. 9 for i? = 5 and R = 7. In the latter 
cases, the true asymptotic behavior is not accessible and is cut off by the finite value of L. 

In addition, a more detailed analysis of the Monte Carlo results for small values of R but large values of L explicitly 
shows that the decay of the bound density deviation behaves as ^ l/\/a; for large x but decays faster than 1/y/x for 
smaller values of x. This crossover behavior is shown in Fig. 10 where the simulation data are compared with the thin 
dashed line corresponding to the decay law ~ l/\/x. Furthermore, the data can be well fitted with a density profile 
as given by (24) if one makes an appropriate choice for the scale parameters Dsc and Xsc- 

The Monte Carlo data shown in Fig. 10 correspond to a tube of radius R = 3 and length L = 6000 with boundary 
densities pb,in = 0.8 and pb.ex = 0.5. In this case, a least-squares fit of the data for a; > 20 leads to the parameter 
values Dsc — 0.81 and x^c — 8.57. The corresponding fitting curve corresponds to the thick dashed line in Fig. 10. 
The fit becomes less reliable close to the left boundary, where the assumption that rjb is small is no longer fulfilled. 
For fixed boundary densities, both fitting parameters Dg^ and Xgc are found to depend on R and to increase with 
increasing R. In addition, for i? > 4, it becomes rather difficult to determine these parameters since one would have 
to simulate rather long systems with L ^ 6000 in order to observe the true asymptotic behavior. 

For the ASEP in one di mensi on, it has been argued that the power-law decay of the density profile as given by 
the asymptotic form « cj \fxjl for large x is characterized by the universal amplitude c = 1/2y^ ~ 0.282. [23,24] 
Inspection of the relations (25), (26), and (28) shows that, in the present situation, c = Dscl2v\,^Xsc where the 
parameters Dsc and Xac depend on the tube radius R. This implies that the amplitude c will also depend on R. Using 
the numerically determined values for D^c and Xgc and the bound state velocity V}, ~ 0.99, we obtain the estimates 
c ~ 0.27 for i? = 0, which corresponds to the 1-dimensional ASEP and should be compared with the exact value 
c — 1/2-^7? ~ 0.282, and c ~ 0.14 for R = 2>. Thus, for small values of R^ the amplitude c is found to decrease with 
increasing R. 
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V. DIFFUSIVE INJECTION AND EXTRACTION OF MOTORS 



Finally, we consider the second type of open boundary conditions corresponding to case (C) in Fig. 2. The length 

of the tube is again denoted by L. This tube is now longer than the filament which has length Lp < L. The left end 
of the tube is located at a; = as before but the left end of the filament is at a; = AL = (L — Lf)/2. Likewise, the 
right end of the filament is at a; = AL + Lp whereas the right end of the tube is at a; = 2AL + Lp = L. Thus, a motor 
particle which enters the tube on the left must diffuse over a distance ^ AL before it can come into contact with the 
filament, and a motor particle which leaves the filament at its right end must also diffuse over a distance ~ AL before 
it can leave the tube. 

At the left and right end of the tube, we now prescribe constant boundary densities as given by 

p{x = Q,y,z) = puh,in and p{x = L + l,y,z) = puh,eK (34) 

for all values of y and z with y'^ + z"^ < . 

As before, the jump probability /3 to make backward steps on the filament is taken to be zero. Likewise, the resting 
probability 7 to make no step at all on the filament is also zero with a possible exception at the 'last' filament site 
with {x,y,z) = (AL + Lf,0, 0). Indeed, in order to define the system in a unique way, we still have to specify the 
probability to make a forward step at this 'last' filament site. Two possible choices appear rather natural: (i) Active 
unbinding from the 'last' filament site, i.e., the motor particle attempts to step forward with probability a and makes 
a step if the adjacent nonfilament site is unoccupied. In this case, the forward step at the last filament site is governed 
by the same probabilities as all other forward steps along the filament; and (ii) Thermal unbinding in which the motor 
particle unbinds with probability e/6 both in the forward direction and in the four orthogonal directions. In the latter 
case, one has to choose the resting probability 7 to be nonzero and to be given by 7 = 1 — 5e/6. 

The choice (i) is suggested by the results of recent experiments on microtubules and kinesin motors [25] which 
indicate that these motors unbind quickly at the filament ends. In the following subsections, we will first consider 
this choice (i) corresponding to active unbinding from the last filament site. In the last subsection VD, we will show 
that the choice (ii) leads to rather similar behavior. 



A. Diffusive bottlenecks 



In order to understand the behavior found for boundary condition (C), it is instructive to partition the tube into 
three compartments which are defined as follows: (i) A left compartment with 1 < a; < AL where transport is purely 
diffusive; (ii) A middle compartment with AL < x < AL + Lp where all directed (or active) transport occurs; and 
(iii) A right compartment with AL + Lp < x < L where the transport is again purely diffusive. 

For a stationary state, the total current through the tube must be constant and, thus, must be the same in all three 
compartments. The current through the middle compartment is given by the bound current j'b = "^hP^i^ ~ Ph) ■ Thus, 
the diffusive currents Jdif,L and Jdif,R in the left and right tube segment must be equal and must satisfy the simple 
relation 

^dif,L = Jdif.R = VbPb(l - Ph) ■ (35) 

The relation as given by (35) is easily checked in the simulations since the density profile is found to be approximately 
linear in the left and in the right compartments provided AL is sufficiently large. For such a linear density profile in 
the right compartment, the diffusive current Jdif,R can be estimated as 

T /^ ^AT \r) |Pub(a; = AL + Lp) - p^b,ex| 

Jdif.R - (1 + -'Vchji'ub • (36) 

Since the maximal density difference is one (in the units used here), the maximal diffusive current behaves as 

max(Jdif,R)~(l + Areh)-Dub^~i?'i?ub^ . (37) 

This shows that the maximal diffusive current depends on the tube radius i?, the lateral size AL of the diffusive 
segments, and the diffusion coefficient D^\y of the unbound motors. Since these parameters can be chosen independently 
from the bound motor velocity Vb, the diffusive current Jdif,L can be made smaller than the maximal bound current 
Vb/4 on the filament. In the latter case, the diffusive compartments act as diffusive bottlenecks and the maximal 
current phase characterized by the current v^/A is expected to be absent from the phase diagram. This expectation 
is indeed confirmed by the simulations as discussed next. 
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B. Phcise diagram without maximal current phcise 

One geometry, for which no maximal current phase has been observed in the simulations, is provided by a tube 

with radius i? = 5, length L = 600, and filament length Lp = 590. The corresponding phase diagram as determined 
by the Monte Carlo simulations is shown in Fig. 11. The largest part of the phase diagram is covered by the high 
density phase; in addition, a low density phase is found for small values of the boundary densities Pub.in and Pub.ex- 

The transition line displayed in Fig. 11 has been determined from the hmctional dependence of the bound density 
in the bulk, pjj, on the left boundary density Pub.in as shown in Fig. 12. Inspection of this latter figure shows that 
the bound bulk density p° jumps at a certain value of the left boundary density Pub,in- Since L is finite, this jump 
occurs over a small but finite interval of Pub.in- Thus, the jump Ccin be cliarcLCtcrizcd by two /?ub,in~ 

values, say p^^ .^^ 

and p^jj jjj, which represent the left and the right 'corner' of the numerically determined jump. Both pub,in-values 
have been included in the phase diagram of Fig. 11. 

For p^^ < Pub.in < Pub in' simulations do not reach a stationary state within two days of computation. 
Simulations also become very slow in the high density phase especially when the overall motor concentration gets 
so large that the bound density is close to one and the unbound density is no longer small compared to the bound 
density. 



C. Presence of the maximal current phase 

In order to estimate the set of parameters, for which the phase diagram exhibits a maximal current phase, we return 
to the estimate (36) for the diffusive current Jdif ,R and make the simplifying assumption that the bound and unbound 
densities in the middle compartment of the tube are essentially constant. Thus, we replace the true unbound density 

puh{x = AL + Lp) at the right end of the filament by its bulk value p'^^. 



The bulk density p'^^ of the unbound motors is related to the bulk density p° of the bound motors via the radial 



equilibrium relation (3). For the maximal current phase, one has p^ = 1/2 and (3) leads to p^j^ = i^^/!^' ~ e/n^.d 



for small e. In this way, we arrive at the estimate 

Jdif.R ^ (1 + iVch)-Dub — ^ for p„b,ex = . (38) 

The maximal current phase should be present in the phase diagram if the diffusive current Jdif.R exceeds the 
maximal current v^/A on the filament. It then follows from (38) that the maximal current phase should be present 
for Pub, ex = provided the tube radius R satisfies 

i?^>i?J.^^ (39) 
47r i^ube 

where 1 + A'ch ~ tt-R^ has been used. Using the same line of reasoning, a second, less restrictive condition can be 
obtained from an estimate for the diffusive current within the left compartment. 

The threshold value R^, for the tube radius as given by (39) has been confirmed by Monte Carlo simulations for 
the jump probabilities a = 1 — 2e/3 with e = 1/100, the sticking probability TTad = 1, the compartment size AL = 5, 
and the boundary densities Pub,in = 0.2 and Pub.ex = 0. The jump probabilities imply the bound motor velocity 
Vb = 1 — 2/300; the diflFusion coefficient D^h of the unbound motors has the value Dub = 1/6 (since the resting 
probability 7 = as mentioned above). When these parameter values are inserted into (39), one obtains the estimate 
i?* ~ 15.4. 

The corresponding Monte Carlo data are displayed in Fig. 13. Inspection of this figure shows that the current does 
indeed attain its maximal value Ub/4 for i? > i?» with i?, ~ 16. The same Fig. 13 also shows the transition from the 
low density to the high density phase which occurs for a tube radius _R»» which satisfies 4 < R^,^, < 5. 

A complete phase diagram for a tube with radius i? = 17 is shown in Fig. 14. Again, most of the phase diagram is 
covered by the high density phase (HD), while a low density phase (LD) is found only for very small values of Pub.in- 
The maximal current phase (MC) is present now but only for very small values of Pub, ex- 
it is interesting to note that similar effects also occur in the purely one-dimensional system if one considers a 
driven system which is bounded by two segments which exhibit only diffusive transport. In this 1-dimensional case, 
quantitative predictions can be made using a mean field approximation as will be discussed elsewhere [26]. 
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D. Active versus thermal unbinding from the 'Icist' filament site 



The Monte Carlo data displayed so far have been obtained for active unbinding of a motor particle which is bound 

to the 'last' filament site. As mentioned above, another possibility is that the motor particle gets stuck at the 'last' 
filament site and unbinds only by thermal excitations, i.e., with unbinding probability e/6. For these two different 
unbinding mechanisms, one will, in general, obtain different density profiles. However, this diflFerence is not dramatic 
as one can see from Fig. 15 which exhibits density profiles for both cases. Although the probability for a forward step 
at the 'last' filament site differs by two orders of magnitude for the two cases, the bulk density exhibits a relatively 
small difference. The bound density, on the other hand increases and decreases close to the right end of the filament 
for thermal and active unbinding from the 'last' site, respectively. 



Let us summarize the main results and add a few comments. We have studied a lattice model for the motion of 
many molecular motors in an open tube which contains a single filament. When bound to the filament, the motor 
particles undergo an asymmetric simple exclusion process (ASEP). In addition, motors can unbind from the filament 
and then diffuse freely in the tube. As for the ASEP in one dimension, the motor traffic in open tubes can exhibit three 
different phases: high density and low density phases which are characterized by an exponential decay of the density 
deviations from their bulk values and maximal current phases characterized by an algebraic decay. Therefore, the 
molecular motor traffic in open tubes are promising candidates for the experimental observation of boundary-induced 
non-equilibrium phase transitions. 

In general, the location of the transition lines is found to depend on the precise choice of the boundary conditions. 
Apart from periodic boundary conditions, case (A), we studied two different boundary conditions (B) and (C) for 
an open tube. In case (B), the bound and unbound densities are kept fixed at the boundaries and satisfy radial 
equilibrium. For this case, the location of the transition lines is independcnit of the model parameters, and the phase 
diagram of the ASEP in one dimension is recovered. In case (C), the active compartment of the tube is bounded 
by two compartments where the transport is purely diffusive. In this latter case, the phase diagram depends on the 
geometry of the tube and on the transport properties in the bound and imboimd motor states. In many cases, the 
maximal current phase is completely suppressed by the coupling to the diffusive compartments which act as bottle 
necks for the transport. 

The theoretical results described here should be accessible to experiments on cytoskeletal filaments and motors. In 
particular, the motor traffic through open tubes as discussed here provides new opportunities to study the transport 
properties of ASEPs by systematic experiments. 



We thank Theo M. Nieuwenhuizen for comments and discussions during the initial stage of this work. 

APPENDIX A: RADIAL EQUILIBRIUM FOR PERIODIC BOUNDARY CONDITIONS 

In the following, we show that Eq. (3), the condition for radial equilibrium, holds exactly in the case of periodic 
boundary conditions using the quantum Hamiltonian formalism [27]. The exact stationary master equation can be 
written in the form 



where H is the 'quantum Hamiltonian' of the stochastic process and \p) is a vector in a product Hilbert space; each 

lattice site is represented by a two-state system with the orthogonal vectors |1) for an occupied lattice site ("spin 
down") and |0) for a vacancy ("spin up"). Because of translational invariance in the direction parallel to the filament, 
there cannot be any radial currents in our system, and the unbound density is independent of the radial coordinate; 
thus we can restrict the analysis to the case of one unbound channel. We denote by \p)k,h and \p)k,uh the state of 
site k of the bound and unbound channel, respectively. Using the general recipe given in chapter 2 of Ref. [27] we 
construct the 'quantum Hamiltonian' H of our system: 



VI. SUMMARY AND CONCLUSIONS 



ACKNOWLEDGMENTS 



H\p) = 0, 



(Al) 



H = Hi + H2 + H3, 



(A2) 



12 



where Hi represents the dynamics of the asymmetric exchision process in the boimd channel, H2 the symmetric 
exclusion process in the unbound channel and H3 the coupling of the two channels. Each term can be written as a 



sum Hi = J2k 

= ink,hVk+i,h - 4,b«fe+i,b) (A3) 

hf^ = r>ub (nfe,ub^^fe+i,ub - Sfc,ub«fe+i,ub + ^^fe,ubrife+i,ub - Sfc_ub«fe+i,ub) (A4) 

^k^ = \ {nk,hVk,nh - Sfc.bSfc.ub) + ^ (^'fe,b^^fc,ub - Sfc,bSfc,ub) • (A5) 



b is a creation operator for a particle at site fc of the bound channel, and s^^ is the corresponding annihilation 
operator. nk,h is the particle number operator at site fc of the bound channel and Vk,h = 1 — '^fc.b- Operators for the 
unbound channel are defined in an analogous manner. The off-diagonal parts of the operators h), represent hopping 
of particles, while the diagonal parts are determined by conservation of probability, see chapter 2 of Ref. [27]. 

We now show that the product measure 

|P> = (8) ( [(1 - Pb)|0)fe,b + Pb|l)fc,b] ® [(1 - Pub)|0)fc,ub + /Oub|l)fc,ub] ) , (A6) 
fc 

is a stationary state, if Eq. (3) holds, i.e. that the radial equilibrium condition implies H\p) = 0. The product measure 
\p) defines a state where the density is pb at each lattice site of the bound channel and Pub at each site of the unbound 
channel and spatial correlations vanish. 

Hi and H2 do not couple the bound and unbound channel, therefore we can consider them separately and refer 
to the result, that for the symmetric as well as for the asymmetric exclusion process, the product measure \p) is 
stationary in the case of periodic boundary conditions [14]. For a proof using the quantum Hamiltonian formalism, 
see chapter 7.1.2 of Ref. [27]: In the case of the ASEP in one dimension, for example, one can easily check, that 

^^ii'\p) — "b('^fe,b — nk+i,h)\p), therefore the summation gives zero for periodic boundary conditions. Hence Hi\p) = 
and ff2|p) = 0.' 

Concerning H^, it is sufficient to consider a single site fc in both channels. Doing some steps, we find 

("fc.bffe.ub - ■Sfe.b'Sfc.ub) \p}k,h\p)k,uh = Pb(l - Pub)(|l)fc,b|0)fe,ub - |0)fe,b|l)fc,ub) (A7) 

and 

(vk,hnk,nh - Sfe,b*fe,ub) |p)fc,b|p)fe,ub = Pub(l - Pb)(|0)fc,b|l)fc,ub - |l)fc,b|0)fc,ub) (A8) 

for the product measure (A6). Hence from 

h'j^^\p)k,h\p)k,uh = (^Pb(l - Pub) - ^Pub(l - Phi) (|l)fc,b|0)fc,ub - |0)fc,b|l)fc,ub) = 0, (A9) 

we obtain Eq. (3), which is the condition for radial equilibrium. 



APPENDIX B: CONTINUUM TWO-STATE MEAN FIELD EQUATIONS FOR OPEN BOUNDARIES 

In this appendix, we give the details of the continuum mean field approximation for the two-state model introduced 
in section IV B. We insert the decomposition pb = Pb + '^b and pub = Pub + ^ub as introduced in (15) into the mean 
field equations (11) and (12) and expand these equations up to second order in the density deviations and Tjub- As 
a result, we obtain 

d d 

-Db^??b + A^ch-Dub^??ub = ^^b(l - 2pb)r/b - Ub??b (Bl) 



and 



d d d"^ 

Vh{l - 2pl)—r]b - ^VhVh-g^ilh - = A?7ub - 5% + (e - 'S-ad)??ub??b, (B2) 
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with 



A = n^{l-pl) + -epl and S = e(l - + f^adP^b- (B3) 

Note that radial equiUbrium for the densities p\, and Pub would imply that the right hand side of (B2) vanishes. 

Let us first consider the case ^ 1/2 which applies to the high and low density phases. In this case, we can neglect 
the second order terms and obtain two equations which are linear in the density deviations. The second equation 
(B2) can then be solved for r^ub which leads to 

, v^{l - 2pl) d £>b 52 



with 

B _ g(l-p;^J+#adP°b 

^-A %^{l-pl) + ~epl ■ 
When this expression is inserted into (Bl), we obtain the first order relation 



(B5) 



for ?7b. We now make the Ansatz r^b ^ cxp(.T/^) which leads to the cubic equation (16) for the decay length ^. 

Now we consider the maximal current phase, i.e. the case = 1/2. In this case the linear terms are zero. 

Furthermore, we neglect terms of order rjb ^ and as can be justified a posteriori since % is found to decay as 

an inverse power of x. Thus, up to leading order, we can ignore the left hand side of (B2), and the right hand side 
of (B2) vanishes. This implies radial equilibrium for the asymptotic decay to the homogeneous solution. Up to this 
order, we find 

BVb S(e-7fad) 2 m'7^ 



for small rih and therefore 



^r,ub = <?^r;b + 0(r?b^). (B8) 



If this latter expression is inserted into (Bl), we obtain 

' dx 

which leads to 



d 

{Db + 5^chl)ub) ^Vb^-Vbvl (B9) 



77b ~ tor large x . (BIO) 

VhX 

For A^ch = 0) this asymptotic behavior is identical to the one found from the mean field approximation for the ASEP 
in one dimension [15]. 
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APPENDIX C: LIST OF SYMBOLS 



a abbreviation for Dgcf^VbXsc used in subsection IV B 2 

(A) periodic boundary conditions, see Fig. 2 

A abbreviation used in Appendix B, defined in Eq. (B3) 

a probability for a forward step of tlie bound motor 

b abbreviation for 1 + .g iVch-Dub / Dgc used in subsection IV B 2 

(B) open tube, radial equilibrium at the boundaries, see Fig. 2 
B abbreviation used in Appendix B, defined in Eq. (B3) 

f3 probability for a backward step of the bound motor 

(C) open tube, diffusive injection and extraction of motors, see Fig. 2 
D diffusion coefficient for the 1- dimensional ASEP 

£>b diffusion coefficient of the bound motor particle 

AL distance between filament end and tube end for boundary conditions (C) 

Dub diffusion coefficient of the unbound motor particle 

Dsc scale parameter for the scale-dependent diffusion coefficient 

e unbinding probability of a bound motor 

e rescaled unbinding probability 2e/3 

r^b {x) deviation of the bound density from the constant bulk value 

?7ub(a;) deviation of the unbound density from the constant bulk value 

g parameter defined by Eq. (20) 

7 resting probability of a bound motor 

J (global) c;m-ront through the tube 

jb local bound current 

jub local unbound current 

•^dif.L diffusive current in the left diffusive compartment for boundary conditions (C) 

•^dif.R diffusive current in the right diffusive compartment for boundary conditions (C) 

L length of the tube 

£ lattice constant on filament 

Lf length of the filament 

Nch number of unbound channels 

TVmo number of motors in the tube which is fixed for periodic boundary conditions (A) 

(p cross section of the tube 

TTad sticking probability for a motor hopping to the filament 

TTad rescaled sticking probability 27rad/3 

r radial spatial coordinate 

R radius of the tube 

i?* minimal tube radius for which a maximal current phase is found in case (C) 

p{x) density in the case of the 1-dimensional ASEP 

Ph{x) density of motors bound to the filament 

Pb,ex right boundary density on the filament 

Pb.in left boundary density on the filament 
constant bulk density on the filament 

Pex right boundary density for the 1-dimensional ASEP 

Pin left boundary density for the 1-dimensional ASEP 

Pmo overall motor concentration 

/9° constant bulk density for the 1-dimensional ASEP 

Puh{x,r) density of unbound motors 

Pub.ex density of unbound motors at the right boundary 

Pub,in density of unbound motors at the left boundary 

Puh unbound constant bulk density 

T basic time unit 

V velocity for the 1-dimensional ASEP 

Vh velocity of bound motor 

Vc velocity of density fluctuations 

Vs domain wall velocity 

X spatial coordinate parallel to the filament 
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Xsc scale parameter for the scale-dependent diffusion coefficient, see Eq. (23) 

a;*,3;*, crossover lengths 

^ decay length for the density deviations 

^0 localization length for the 1-dimensional ASEP 

xq spatial location of the left tube end 

y rescaled and shifted spatial coordinate, y = {x — xo)/xsc 

y* , y** rescaled crossover lengths 
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FIG. 1. Motor particles which can bind and unbind to a filament (dark rod) within a cylindrical tube. Single motors which 
are bound to the filament (and are not sterically hindered by other motors) move with velocity Vb to the right. Unbound 
motors diffuse with diffusion coefficient D^h in the surrounding liquid. 





(C) 




FIG. 2. Different types of boundary conditions: (A) Periodic boundary conditions which is similar to a closed torus geometry; 

(B) Open tube with boundaries satisfying radial equilibrium: the bound and unbound motor densities are fixed at the two 
boundaries and satisfy radial detailed balance at each boundary; and (C) Open tube with diffusive injection and extraction of 

motors. In all cases, the tube has total length L; in cases (A) and (B), the filament has the same length as the tube; in case 

(C) , the filament has length Lp < L and there are two boundary compartments of linear size AL. 



18 



J/v. 




N /L 
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FIG. 3. Reduced current J/ub through the tube with periodic boundary conditions as a function of the reduced particle 
number TVmo/i. The hue is calculated from Eq. (7), the Monte Carlo data are obtained for a tube of length L — 200 and 
radius R = 25 corresponding to channel number TVch ~ 1940. The random walk probabilities are 13 — 0, j — 99/100, e — 10"'', 
a = 1 — 7 — 2e/3, and vTad = 1. 
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Pb,in 

FIG. 4. Phase diagram for motor traffic in open tubes with boundary condition (B) as a function of the boundary densities 
Pb.in and pb,cx at the left and right end of the filament, respectively. There are three phases distinguished by the bulk value of 
the bound density pJJ: a low density phase (LD), a high density phase (HD), and a maximal current phase (MC). This phase 
diagram is identical to the phase diagram of the 1-dimensional ASEP as explained in the text. 
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FIG. 5. Bound and unbound density profiles pb and p^b as a function of the spatial coordinate x for the low density phase 
with boundary conditions (B) and boundary densities pb.in ~ 0.3 and pb.ex ~ 0.6. The unbound density has been averaged 
over the tube cross section and multiplied by a scale factor of 10. All data points are averages over six lattice sites in direction 
parallel to the filament. The tube has length L = 600 and radius R = 10. The random walk probabilities are /3 = 7 = 0, 
e = 10~^, a = 1 — 2e/3, and TTad = 1. 
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FIG. 6. Bound and unbound density profiles pb and Pub as functions of the spatial coordinate x for the high density phase; 
boundary conditions (B) with pb,in = 0.8 and pb,ex = 0.6. The unbound density has been averaged over the tube cross section 
and multiplied by a scale factor of 10. The geometry and the random walk probabilities are the same as in Fig. 5. 
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FIG. 7. Localization length ^ as a function of the tube radius R for boundary conditions (B) with pb.in = 0.38 and pb.ex = 0.6. 
Circles arc simulation data, obtained for a tube with length L = 600 and with random walk probabilities as in Fig. 5; lines are 
the corresponding results of mean field calculations (dotted: discrete two-state approximation, dashed: continuous two-state 
approximation, solid: full diffusion equation, see text). The two crosses and the filled circle at i? = represent the mean field 
results and the exact result for the 1-dimensional ASEP, respectively. 
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FIG. 8. Bound and unbound density profiles pb and pub as functions of the spatial coordinate x for the maximal current 
phase; boundary conditions (B) with pb,in = 0.6 and pb,ex = 0.1. The unbound density has been averaged over the tube cross 
section and multiplied by a scale factor of 10. The parameters are the same as in Fig. 5. 
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FIG. 9. Deviation ijb of the bound density from its constant value fax from the boundaries as a function of the spatial 
coordinate x for the 1-dimensional ASEP (solid line) and for motor traffic in tubes with radius -R = 3, 5, 7 (dashed lines) and 
boundary conditions (B) with pb.in = 0.8 and pb.ex = 0.5. The curves for the 1 dimensional ASEP and for the tube radii 
R = 3,5 have been multiplied by scale factors 0.1, 0.2, and 0.5, respectively. The thin dotted lines correspond to the decay 
laws ~ l/a and ~ l/^/x, respectively. The tube length is L = 2000 and the random walk probabilities as in Fig. 5. 
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FIG. 10. Deviation jjb of the bound density from its constant value far from the boundaries as a function of the spatial 
coordinate x for a tube with radius R = 3 and length L = 6000; boundary conditions (B) with Pb,in = 0.8 and pb.cx = 0.5. The 
thin solid line is obtained from the simulation, the thick dsished line is a fit to Eq. (24) with Dsc — 0.81 and Xsc — 8.57. The 
dotted line indicates the power law ~ l/-\/x. Note that the simulation data decay faster than ~ l/\/S for lO'^ ^ x ^ 10^. The 
random walk probabilities are the same as in Fig. 5. 
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FIG. 11. Phase diagram for tubes with boundary conditions (C) as a function of the left and right boundary densities pub, in 
and Pub,ex for a set of parameters where no maximal current phase occurs. The tube has length L = 600 and radius R = 5, 
the filament length is Lf = 590 and the distance between the filament ends and the tube ends is AL = 5. The random walk 
probabilities are the same as in Fig. 5. Note the different scales of the axes. 




0.2 0.4 0.6 0.8 1 
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FIG. 12. Bound density function of the left boundary density Pub,in for boundary conditions (C) and different values 

of the right boundary density Pub.ex. The parameters are the same as in Fig. 11. The discontinuity in the functional dependence 
of Pb on Pub, in corresponds to the transition from the low to the high density phase, compare Fig. 11. 
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FIG. 13. Current J/vh (diamonds) and bound density (circles) as a function of the tube radius R for boundary conditions 
(C) with Pub, in = 0.2, Pub, ex = 0. The tube length L — 600, the filament length Lp — 590, and the distance between the 
filament ends and the tube ends is AL — 5. The random walk probabilities are the same as in Fig. 5. 
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FIG. 14. Phase diagram for a tube of radius R — 17 with boundary conditions (C) as a function of the left and right boundary 
densities Pub.in and pub,cx. The other parameters are the same as in Fig. 11. In the inset, which shows the complete phase 
diagram, the maximal current phase can be hardly distinguished from the line pub,Gx = 0. 



24 






. bound 








unbound 

( coccoscccccccaxxxxx 





0.8 



0.6 



0.4 



0.2 



50 550 



560 



570 



580 



590 



600 



X 



FIG. 15. Bound and unbound density profiles pb and pub as functions of the spatial coordinate x for boundary conditions 
(C) with Pub, in = 0.2 and pub,ox ~ 0. Two different unbinding processes of the motor particles at the 'last' filament site are 
compared: (Filled circles) Active unbinding for which the motor makes one final step in the forward direction; and (Open 
Circles) Thermal unbinding with probability e/6 in the forward direction. The geometric parameters are L = 600, Lp ~ 590, 
AL = 5, and R = 5; the random walk probabilities are as in Fig. 5. 
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